Self-limited oxide formation in Ni(lll) oxidation 
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The oxidation of the Ni(lll) surface is studied experimentally with low energy electron microscopy 
and theoretically by calculating the electron reflectivity for realistic models of the NiO/Ni(lll) 
surface with an ab initio scattering theory. Oxygen exposure at 300 K under ultrahigh-vacuum 
conditions leads to the formation of a continuous NiO(lll)-like film consisting of nanosized domains. 
At 750 K, we observe the formation of a nano-heterogeneous film composed primarily of NiO(lll) 
surface oxide nuclei, which exhibit virtually the same energy-dependent reflectivity as in the case of 
300 K and which are separated by oxygen-free Ni(lll) terraces. The scattering theory explains the 
observed normal incidence reflectivity R(E) of both the clean and the oxidized Ni(lll) surface. At 
low energies R(E) of the oxidized surface is determined by a forbidden gap in the ky = projected 
energy spectrum of the bulk NiO crystal. However, for both low and high temperature oxidation 
a rapid decrease of the reflectivity in approaching zero kinetic energy is experimentally observed. 
This feature is shown to characterize the thickness of the oxide layer, suggesting an average oxide 
thickness of two NiO layers. 

PACS numbers: 68.37.Nq, 81.65.Mq, 71.15.Ap 
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I. INTRODUCTION 

Surface oxidation is an almost ubiquitous phenomenon, 
and it is generally associated with profound changes in 
geometrical structure and materials properties. The late 
transition metals (TMs) have received persistent atten- 
tioni owing to their tremendous importance in a vari- 
ety of hetcrogencously catalyzed chemical reactions^ The 
general oxidation mechanism may be qualitatively quite 
successfully described by the Cabrcra-Mott model^ devel- 
oped more than 60 years ago, from which many common 
oxidation phenomena, e.g., the temperature-dependent 
thickening of oxide films until saturation, can be ratio- 
nalized^ However, the different elemental properties of 
the TMs, arising from the varying d-band occupation 
through the elemental series, give rise to a huge diver- 
sity in oxidation pathways and here, using the oxida- 
tion of Ni as an example, we present a case where even 
this fundamental concept of oxide saturation thickness 
increasing with temperature needs to be rephrased. 

A prominent example is the interaction of molecular 
oxygen with Ni(lll), which is frequently viewed as a 
model system for dissociative adsorption of molecular 
oxygen on a TM surface. Numerous studies have al- 



ready targeted the initial adsorption of oxygen onto clean 
Ni(lll), which induces a (2x2) surface reconstruction.^ 

While the structural model of this low-coverage phase, 
which is established upon low O2 exposure, is well ac- 
cepted, there is no general consensus regarding the struc- 
ture of the nickel oxide layers that form with prolonged 
oxygen dose. An especially interesting case is the range 
of sample temperatures above 600 K during oxidation, 
which has so far remained mostly unexplored. Below 
600 K, different scenarios have been proposed depend- 
ing on temperature, involving the evolution of few-layer- 
thick, bulk-like NiO filmsZr— or the formation of individ- 
ual oxide grain o 11 ' 12 of thicknesses exceeding 10 nm^ 

Due to this inherent complexity of surface oxidation at 
its early stages, deeper insight into the underlying physi- 
cal and chemical processes may be gained from informa- 
tion gathered in situ, i.e., acquired in real time during 
reaction in an oxidative environment. Ideally, the ex- 
perimental tools of choice should enable a simultaneous 
surface-sensitive characterization of geometrical and con- 
comitant electronic structure, establishing a link between 
the two properties. 

A promising approach to real-time studies of dynamic 
surface processes in TM oxidationi^r— is to employ scat- 
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tcring of low-energy electrons in full- field microscopy, also 
known as low-energy electron microscopy (LEEM) i 19 ' 20 
This method is intrinsically sensitive to the near-surface 
region owing to the strong interaction of slow electrons 
with condensed matter, and it provides microscopic in- 
sight on a length scale of a few nanometers. Infor- 
mation on the surface crystal structure is accessible by 
the related technique of low-energy electron diffraction 
(LEED), which allows to determine the dimension of the 
surface unit cell and the point group of the lattice. In ad- 
dition, for a given reflected beam the dependence of the 
reflected current / on the acceleration voltage V of the in- 
cident electrons, the /(V) curve, contains information on 
the structure of the unit cell. The atomic arrangement on 
the surface can be inferred from the multiple-scattering 
analysis of the I(V) curves for several reflected beams 
over a wide energy rangei 21 > 22 This technique, commonly 
referred to as intensity- voltage, /(U)-LEED, has already 
enabled the determination of many surface structures^ 
The /(V) curves also reflect the bulk electronic structure, 
in particular, the Bragg gaps and critical points, where 
the band structure strongly deviates from free-electron- 
like behaviori 24 i 25 This effect is especially pronounced at 
low energies, where the inelastic scattering is rather weak. 

Experimentally, the much higher intensity of the re- 
flected electron current at very low energies in the range 
of a few eV, typical of LEEM measurements, makes it 
very attractive to record the I(V) dependence in a spa- 
tially and time-resolved manner by sampling the elec- 
tron energy in imaging mode. This so-called /(V)-LEEM 
technique facilitates the measurement of the individual 
I(V) dependence of nanosized surface phases and has 
already been applied to processes in epitaxial growth,— 
surface chemical reactions^ and oxidation catalysisi 16 i 17 

However, despite the achievements of the well- 
established techniques in the interpretation of LEED at 
the energies of a few tens of eV, their application at very 
low energies is not straightforward: most of the imple- 
mentations of the multiple scattering theory rely on a 
rapid decay of the electron wave into the solid, while at 
the energies of a few eVs the inelastic scattering is weak, 
and the penetration depth is very large. In addition, 
low-energy electrons are rather sensitive to details of the 
crystal potential^ so for a fully conclusive comparison 
with the experiment the LEED calculation should rely 
on the self-consistent potential - in the sense of density 
functional theory (DFT) - both in the bulk and at the 
surface. 

In the present work we present a combined experimen- 
tal and theoretical study of the energy-dependent con- 
duction properties of oxygen overlayers on Ni(lll) with 
the aim to relate the observed I(V) curves to the elec- 
tronic properties of the surface and, eventually, to its 
atomic structure. We apply an ab initio scattering theory 
based on a full-potential augmcnted-plane- waves (APW) 
formalism to a model of the oxidized Ni(lll) surface. 
A good agreement between the calculated and the mea- 
sured LEED spectra allows us to interpret the spectral 



structures and to conclude on the thickness of the sur- 
face oxide layer. In particular, we find a forbidden en- 
ergy gap in the k|| = projected spectrum of the (111) 
surface of the bulk NiO crystal, which causes a rapid in- 
crease of electron reflectivity from the oxidized Ni(lll) 
surface at low kinetic energies. We show that with de- 
creasing thickness of the oxide layer the electron trans- 
mission does not increase uniformly over the gap region, 
but a narrow transmission channel opens at the bottom 
of the gap, which is observed in our /(V) measurements. 

The paper is organized as follows: In Sees. Inland Mil 
we describe the experimental procedures to oxidize the 
pristine Ni(lll) surface and characterize the oxygen-rich 
phases in situ. The theoretical approach is introduced in 
Sec. IIV1 and applied to Ni(lll) and to thin NiO layers 
on Ni(lll) in Sees. fVl and IVI1 respectively. 



II. EXPERIMENTAL DETAILS 

Most experiments were conducted at the National Syn- 
chrotron Light Source (NSLS) at Brookhaven National 
Laboratory (BNL), Upton, NY (USA) using the spec- 
troscopic photocmission and low-energy electron micro- 
scope (LEEM III including hemispherical energy ana- 
lyzer, Elmitec) installed at beamline U5UAj22. Additional 
LEEM measurements were performed in the newly- 
installed LEEM III system (Elmitec, no energy filter) at 
our home institute at the University of Bremen. Trans- 
ferability of the results was cross-checked and asserted by 
performing experiments under nominally identical condi- 
tions. 

A commercially purchased, polished Ni(lll) single 
crystal (Matcck) with a nominal orientation better than 
0.1° was used. After insertion into the UHV chamber 
the sample was cleaned by several cycles of 0.5 keV Ar + 
ion sputtering followed by thermal annealing at 1050 K. 
In addition, short flashes to 1300 K were found to im- 
prove the smoothness of the surface on the nano- to 
sub-micrometer scale. Surface cleanliness following this 
recipe was already asserted in an earlier study^ Sam- 
ple temperatures are given based on the reading of a 
W/Re thermocouple permanently attached to the sample 
support. In the oxidation experiments, research-grade 
(99.998%) oxygen (Matheson Tri-Gas Co.) was dosed 
from a high-precision leak valve in back-filling mode. The 
BNL-LEEM and the Bremen-LEEM systems exhibited 
base pressures of 2 x 10~ 10 Torr and 7 x 10~ u Torr, 
respectively. 

The kinetic energy of the incident electrons in all I(V) 
curves presented in this article is referenced to the onset 
of the mirror electron mode, i. e., with respect to the vac- 
uum level. Experimentally, this onset has been measured 
by determining the inflection point of the abrupt edge in 
sample reflectivity observed when increasing the electron 
energy to a few eV. 
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III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Ni oxidation and NiO structure: status quo 

We begin by briefly reviewing the present status quo 
for oxidation of Ni(lll). Depending on oxygen dose, 
Ni(lll) has been found to undergo a (2x2) surface re- 
construction for coverages up to 0.25 monolayer (ML) 
[for Ni(lll), 1 ML corresponds to an atomic density of 
3.72 x 10 15 cm -2 ] followed by nucleation of an epitaxial 
oxide layer for larger oxygen exposures^— At room tem- 
perature, oxidation has been reported^i^ to start upon 
an accumulated dose larger than 10 Langmuir (L) (1 L 
amounts to 10 -6 Torr-s). A three-stage model of the 
process has been proposed^^I involving (i) dissociative 
chemisorption, (ii) nucleation of oxide islands exhibiting 
a thickness of a few atomic layers, and (iii) slow thicken- 
ing of the oxide film. 

Depending on sample temperature during oxida- 
tion, different NiO orientations have been identified. 
Below 470 K, only NiO(lll) formation has been 
observed ) 5 ' 9 ' 10 ! 32 while NiO(OOl) islands have been found 
for more elevated temperatures^ ' 32 ' 33 These findings 
have been rationalized^ 3 - based on the thermodynamic 
instability of the unreconstructed NiO(lll)-(lxl) sur- 
face, which, because of its rocksalt structure and within 
a picture of purely ionic bonding, should carry a di- 
verging dipole moment^ Since the NiO (111) films were 
experimentally determined to be very thin, i.e., only 
a few atomic layers thick, a metastable character of 
the NiO(lll) orientation has been postulated, in agree- 
ment with the growth of the thermodynamically sta- 
ble, unpolar NiO(OOl) face at 500 K, i.e., at sufficiently 
high temperature^ In other studies, the stability of the 
NiO (111)- (lxl) surface has been related to the adsorp- 
tion of hydroxyl specie o 32 ' 35 that prevent the surface from 
undergoing a so-called "octopolar" (2x2) reconstruction 
leading to a compensation of the surface dipole i 36 ' 37 
Since hydroxyl species are easily removed by annealing to 
600 K)2£ this scenario would provide an alternative expla- 
nation for the stability of the (lxl) phase at room tem- 
perature and the preferential growth of NiO(OOl) above 
500 Ki^2 Yet, it has still remained unclear whether the 
driving force for the (2x2) reconstruction at high temper- 
ature would be the same for thin 3 ^ as well as moderately 
thick NiO (111) films — a question raised almost twenty 
years agod 

Quantitative information on actual oxide thickness, 
however, has proven difficult to obtain, and substantial 
variations in oxide "saturation thickness" have been re- 
ported depending on both temperature and oxide phase. 
At 300 K, several studies have concluded on an ox- 
ide thickness of a few monolayers for the NiO(lll) 
layer jLI .! 3 ^ 3 . In the case of NiO(OOl), this saturation 
thickness may extend to more than 10 nm according to 
medium-energy ion scattering (MEIS)±2 after exposure 
to O2 at 500-600 K, which has been noticed to grow at a 



considerably higher rated 

Based on MEISd it has also been postulated that, in 
the elevated oxidation temperature range, the NiO/Ni 
interface is very rough, consisting of irregularly shaped 
NiO grains, and indeed multiple oxygen-enriched grains 
have been seen in photoelectron microscopy with oxygen- 
sparse areas in between*^ However, until the present the 
high-temperature oxidation regime above 500-600 K has 
remained largely unexplored with respect to atomic and 
nanoscale structured 

In the following Sees. IIII Bl and IIII CI we present our 
results for Ni(lll) oxidation at 300 K and 750 K. In the 
former case, we find the formation of a tight patchwork 
of nanosized NiO(lll) domains. The corresponding I(V) 
spectrum exhibits a very strong, characteristic peak at 
about 6 eV. At a sample temperature of 750 K, which is 
significantly higher than the temperature range explored 
so far, we observe essentially the same local I(V) curve 
for the oxidized parts of the Ni(lll) surface concomitant 
with a local (lxl) LEED pattern, corroborating the nu- 
cleation of ultrathin (lll)-oriented oxide domains even 
at these high oxidation temperatures. 



B. Oxidation pathway at room temperature 

Following the cleaning recipe outlined in Sec. [TT1 a 
clean, well-ordered Ni(lll) surface is established exhibit- 
ing only steps and step bunches in LEEM (sec Fig.QJa)). 
Surface crystallinity is asserted by a sharp (lxl) LEED 
pattern (Fig. Efa)). Afterwards, the sample was exposed 
to molecular oxygen (O2) at a background pressure of 
7.5 x 10~ 8 Torr while the modification of the surface was 
monitored by LEEM (Fig. [TJ and LEED (Fig.©. After 
dosing O2 for 13 s, which corresponds to an accumulated 
dose of 1 Langmuir (L), the integral sample reflectiv- 
ity has considerably decreased, as can be deduced from 
Fig.[IJb). This change in intensity is accompanied by the 
advent of a (2x2) reconstruction in LEED (Fig. Hfb)), 
which is expected for a nominal oxygen coverage of about 
0.25 ML£ 

However, in addition to the integral change in reflec- 
tivity, we observe significant local variations of the (00) 
intensity at structural defects (see Fig.QJb)) that do not 
contribute to the periodic part of the LEED pattern. 
While the terraces exhibit a rather low reflectivity at the 
chosen kinetic energy of the incident electrons (8.0 eV), 
the steps appear significantly brighter in contrast. Since 
these white regions eventually fill the entire surface upon 
prolonged exposure (Fig.QJc))^ we may attribute these 
bright areas to the formation of a more oxygen-rich sur- 
face phase that is associated with a distinct LEED pat- 
tern (Fig. ©c)) of apparent six-fold symmetry but in- 
creased lattice parameter as compared to clean Ni(lll). 
Further dosing of O2 does not induce any visible changes, 
neither in the reflectivity nor in the sample morphology. 
Nevertheless, we note an overall grainy appearance of 
the surface, suggesting the presence of structural features 
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FIG. 1. (color online) (a-c) Low-energy electron micrographs 
showing the transformation of the Ni(lll) surface at room 
temperature upon exposure to molecular oxygen for L (a), 
1 L (b), and 265 L (c). (d) I(V) spectra of clean Ni(lll) 
(dotted line) and transformed regions (solid line). 



just below the resolution limit of the microscope. 

In the following sections, we will mainly focus both 
experimentally and theoretically on the structural iden- 
tification of this bright phase. A comparison of I(V) data 
(Fig-Did)) acquired for the clean Ni(lll) surface with the 
bright, homogeneous area readily corroborates the notion 
of a distinct phase. Further insight into the crystallo- 
graphic structure is accessible by the quantitative analy- 
sis of line profiles extracted from LEED patterns recorded 
during thermal annealing from 300 K to 820 K. While 
no LEED spots associated with the Ni(lll) substrate 
are found after preparation at 300 K, which confirms 
the presence of a continuous, oxygen-rich "film" , these 
integral order reflections re-appear at elevated tempera- 
ture and clearly dominate the LEED pattern at 700 K 
(Fig. [UJd)), whereas the spots of the oxide film have al- 
most vanished. Interestingly, additional spots emerge 
at half-order positions, underlining the formation of a 
(2x2)-0 adlayer phase upon film disintegration and con- 
comitant thermal desorption of oxygen atoms. A quanti- 
tative evaluation using the post-anneal pattern as refer- 
ence reveals a lattice parameter of 4.14 ± 0.03 A for the 
film, which is in good agreement with the lattice constant 
of bulk NiO(lll) (a = 4.177 A)^i indicating an almost 
fully relaxed NiO film. Furthermore, from the width of 
the (lxl) ox peaks at 300 K, an average domain size of 
approximately 1.3 nm can be deduced. 

In separate experiments, room temperature grown 
NiO(lll) films were annealed in a step- wise manner un- 
der other identical conditions in LEED as well as LEEM 
mode. With increasing temperature, the oxide peaks be- 
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FIG. 2. (a-d) Low-energy electron diffraction patterns ob- 
tained after exposing the Ni(lll) surface at room temperature 
to an oxygen dose of (a) L, (b) 1 L, (c) 265 L, and (d) af- 
ter subsequent annealing to 700 K. (e) Evolution of the line 
profile extracted from a cut through reciprocal space along 
the direction indicated in (c, d) with annealing temperature. 
Note that the line profiles were not corrected for variations in 
detector efficiency and that they were taken from a different 
data set than images (c-d). 



come slightly sharper while losing intensity (Fig. EJe)). 
Comparable integral intensities of the oxide and sub- 
strate LEED spots are observed at about 620 K. This 
behavior suggests an Ostwald-like ripening process, in 
which smaller patches cither dissolve or coalesce to form 
larger oxide islands, which yet remain too small to be im- 



aged in LEEM mode. These findings confirm the results 
of a scanning tunneling microscopy study that reported 
an increase in domain size by about a factor of two after 
annealing to 700 Finally, at a temperature of 820 K, 
these islands have completely disappeared, and the orig- 
inal LEED pattern of the clean substrate is restored. 

Summarizing, the presented LEED data strongly 
suggest the formation of a continuous NiO(lll)-likc 
film, which is composed of many small domains of 
a few nanometers in diameter. The recorded I(V) 
curve provides more information on the local crystallo- 
graphic structure, which will theoretically be addressed 
in Sec. IVI1 but at this point it may already serve as a 
fingerprint for a NiO(lll)-likc phase in the upcoming 
oxidation studies performed at elevated temperature. 



C. Oxidation at 750 K 

The LEEM time-lapse sequence for oxygen exposure at 
a sample temperature of 750 K is depicted in in Fig. E2 
On the clean surface (Fig. E3! a )) only steps and step 
bunches arc visible. However, after a 9 L dose we note a 
qualitatively different behavior as compared to the pre- 
vious case of oxidation at room temperature. Under the 
present conditions, islands exclusively nucleate at step 
bunches whereas the single steps and flat terraces remain 
unaffected. Apart from any morphological objects on the 
terraces, we do not find a large-scale change in reflected 
intensity as we observed for oxidation at room tempera- 
ture. This qualitatively different result is confirmed by 
comparing the /(V) fingerprint of the terraces after an 
exposure of about 1000 L (Fig. HJc)), e.g., extracted at 
point "A" in Fig. HJa), with the I(V) reference spectra 
for the clean Ni(lll) surface (Fig.QJd)). Since the curves 
are virtually identical, this result clearly shows that only 
a negligible amount of chemisorbed oxygen is found on 
the terraces at this high temperature. Given the compar- 
atively large oxygen doses needed to achieve substantial 
surface coverage, the integral dissociative sticking coeffi- 
cient for molecular oxygen has to be substantially lower 
than at room temperature, in agreement with previous 
studies that targeted Ni(lll) oxidation kinetics^ 

The few existing oxide nuclei, however, grow upon pro- 
longed exposure and gradually spread over the remaining 
Ni(lll) terraces, while the crossing of step bunches is not 
observed. We also note that essentially no additional, iso- 
lated nuclei appear on a 10 (im scale (Fig. (3jb-e)), which 
is indicative of a significantly enhanced diffusion length of 
the oxygen species. Interestingly, the reacted areas are 
by no means single- crystalline domains of nickel oxide, 
but instead are composed of a variety of small grains ex- 
hibiting different contrast and irregular shapes. From the 
time- lapse sequence (Fig. |3fb-f)), we conclude that the 
edges of the already oxygen-rich areas serve as effective 
nucleation centers for subsequently grown oxide islands, a 
phenomenon that has already been observed in oxidation 
studies of other TM surfaces, e.g., Ru(0001)J£ Because 
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FIG. 3. Low-energy electron micrographs showing the trans- 
formation of the Ni(lll) surface upon oxygen exposure at a 
sample temperature of 750 K. (a) L, (b) 9 L, (c) 40 L, (d) 
215 L, (e) 420 L, and (f) 970 L. 



the lateral growth rate remains steady but slow through- 
out the remainder of the integral oxygen dose of about 
1000 L (Fig.^f)), a considerable fraction of surface area 
remains unreacted to oxygen, in accordance with the in- 
tensity variations found in oxygen concentration maps at 
temperatures of 573 K and 673 K^i 

In Fig. B|a), three main contrast levels are observed, 
one for the oxygen-free Ni(lll) terraces and two in the 
oxide-covered areas. Hence, it would appear that mainly 
two types of oxide domains have nucleated, labeled "B" 
and "C" in the figure. The magnified view shown in 
Fig. EJb), which is a blow-up of the center region indi- 
cated in Fig. EH a), illustrates that the transformed ar- 
eas are composed of a tight network of nanosized oxide 
patches exhibiting different contrast, but whose struc- 
tural properties are nevertheless amenable to /(V) anal- 
ysis. A comparison of the individual I(V) spectra of 
these oxygen-rich phases, using phases "B" and "C" in 
an exemplary fashion, reveals that they are closely re- 
lated since they all show a very similar characteristic res- 
onance at about 6 cV (Fig. EH ))- Moreover, the I(V) 
curve for "B" is virtually identical to the one observed 
for oxidation at room temperature (Fig. EUd)). Thus, 
we may already speculate at this point that all types of 
patches are some form of NiO(lll), but apparently differ 
with respect to certain structural details. 

Further insight is gained from selected-area LEED 
patterns acquired for all regions individually. While 
the local diffraction pattern for "A" (Fig. E2c)) exclu- 
sively exhibits the expected (lxl) spots of the clean 
Ni(lll) surface, the micro-LEED pattern for region "B" 
(Fig. Eflcl)) shows a well-defined, three-fold (lxl) peri- 
odicity, whose surface unit cell is rotated by 30° with 
respect to the substrate lattice. A quantitative analy- 
sis of the peak positions indicates a lattice mismatch of 
19%, in very good agreement with the presence of a com- 
pletely relaxed NiO(lll) film that shares a different reg- 
istry with the substrate as compared to our observations 
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at room temperature. Additionally, for relatively low 
electron energies we also noticed weak, very broad, and 
threefold-symmetric facet spots (not shown), which will 
be subjected to temperature-dependent investigations us- 
ing micro-LEED and high-resolution LEED4 3 - Interest- 
ingly, we observed qualitatively similar diffraction pat- 
terns for both region types "B" and "C" that only dif- 
fered by the relative intensities of the NiO(lll) spots and 
the broad facet streaks. These findings suggest the pres- 
ence of a tight, complicated mosaic structure consisting of 
well-aligned and, probably, considerably tilted NiO(lll) 
domains for both types of regions, albeit exhibiting dif- 
ferent areal ratios of the untilted versus the tilted regions. 
In this structural model, the measured shape of the /(V) 
curve for region "C" would then represent the sum of the 
individual I(V) curves of minuscule, untilted NiO(lll) 
regions as well as of tilted NiO(lll) regions, with the 
latter barely contributing to the integral intensity of the 
specular (00) beam. 

Previous studie s 32 ' 33 argue that for growth tempera- 
tures of about 500 K, only the (001) orientation of NiO 
should be present. Using our in-situ technique, we indeed 
observe NiO(OOl) formation in an intermediate temper- 
ature range of about 500 K (Ref . [43|) , but based on the 
micro-LEED evidence in connection with the strong sim- 
ilarity of the I(V) curves presented here we clearly show 
that this is not the case for even higher temperatures of 
about 750 K, where NiO(lll) is again prevalent. These 
results are in agreement with previous findings of a higher 
thermal stability of the NiO (111) films as compared to 
the NiO(OOl) in annealing experiments, in which the lat- 
ter was shown to decompose upon annealing to 550 K42 

It should be stressed, however, that we did not observe 
any signs of a (2x2) reconstruction of the oxide domains, 
clearly demonstrating that a NiO(lll)-(lx 1) structure 
that is rotated by 30° with respect to the Ni(lll) sub- 
strate lattice is still stable at a temperature of 750 K. 
Evidently, the change in registry adds to the stability of 
these NiO(lll) domains at temperatures close to decom- 
position of the room temperature grown oxide. 



IV. COMPUTATIONAL METHODOLOGY 

There exist two major approaches to an ab initio treat- 
ment of electron scattering by the surface. The more tra- 
ditional multiple-scattering Green's function metho d 21 ' 22 
employs a representation of the crystal by a finite num- 
ber of atomic monolayers. It has an advantage of avoid- 
ing the calculation of the partial Bloch waves inside the 
crystal and proceeds immediately to the scattering so- 
lution.— This method is especially efficient within the 
muffin-tin approximation (MTA) for the crystal potential 
— although full-potential Green's function methods have 
existed for many years^ the MTA is still widely used in 
the theory of LEE] ). "' In the present work we use the 
alternative Bloch wave approach: the LEED function in- 
side the crystal is sought as a linear combination of the 



a) 








M 







c) • 
| 




• 


47.5 eV 



d) 

« 


| 

• 




27.5 eV 




10 15 20 25 30 
Electron Energy (eV) 



FIG. 4. (color online) (a) Low-energy electron micrograph 
recorded after Ni(lll) oxidation at 750 K. (b) Close-up view 
of the center part of the LEEM image shown in (a) with ad- 
justed color map. (c-d) Selected-area LEED patterns acquired 
from regions "A" (c) and "B" (d). (e) Local I(V) curves ex- 
tracted at points "A", "B", and "C" (shown in (a)) after 
Ni(lll) oxidation at 750 K. 



partial waves, which facilitates the interpretation of the 
LEED spectra in terms of the band structure of the sub- 
strate. 

To calculate the reflected intensities R(E) the LEED 
wave function is obtained as a solution of the Schrodinger 
equation for a semi-infinite crystal, see Fig. [5] The scat- 
tering wave function $ is defined by its energy E and the 
incidence direction of the electron beam. In the plane 
parallel to the crystal surface $ obeys the Bloch theorem 
and is characterized by the 2D Bloch vector k||. In the 
vacuum, far from the crystal surface, it is a superposition 
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FIG. 5. (color online) Lower panel: potential profile at the 
clean Ni(lll) surface (thin red line) and with an oxygen over- 
layer (thick black line). Upper panel: LEED wave function 
(the surface Fourier component Gii =0) for E — Ef = 15 eV 
for the clean Ni(lll) (thin lines) and for the overlayer (thick 
lines). The wave functions in the scattering region are shown 
by dashed lines. 

of the incident plane wave and reflected (propagating and 
evanescent) plane waves. Deep in the crystal the poten- 
tial is periodic, and in the absence of inelastic scattering 
(electron absorption) the partial waves are propagating 
(real k±) and evanescent (complex fcj_) Bloch waves that 
comprise the complex band structure of the semi-infinite 
crystal^ Then, the scattering problem consists in find- 
ing the coefficients of the partial waves. 

The evanescent waves carry zero current, and the elec- 
tron beam is completely reflected whenever E falls in 
a k| | -projected band gap. To take into account inelas- 
tic processes, which reduce the reflectivity, an imaginary 
term, the optical potential — iV\, is added to the potential 
in the crystal half-space. The energy E is kept real, so 
the absorbing potential leads to a spatial damping of the 
wave functions, i.e., Bloch vectors of originally propa- 
gating waves acquire an imaginary part. Hence, electron 
absorption is allowed for, and even in the energy gaps of 
the bulk band structure there is no complete reflection 
anymore^ On average, V\ increases with energy, and an 
approximate dependence Vi (E) can be inferred from the 
curvature of the measured I(V) curvei 28 i 49 In the present 
case we used a linear Vi(E) function chosen so as to ap- 
proximately reproduce whenever possible the sharpness 
of the R{E) peaks. An example for Ni(lll) is shown in 
the inset of Fig. H^d). 

The computation starts with constructing self- 
consistent potentials in the bulk crystal and at the sur- 
face within the local density approximation (LDA) of 
the DFT. The band structure is calculated with the ex- 
tended linear augmented plane wave method (EL APW) , 
using the full-potential augmented Fourier components 
technique^ The potential at the surface is determined 
by a repeated-slab calculation. The slabs comprise nine 
atomic layers and are separated by a vacuum region of 
16 a.u., see Fig. [SJ No structure optimization is per- 



formed: all atoms occupy the positions of the ideal Ni or 
NiO lattice, see Fig. H 

The partial Bloch waves tpk± arc obtained as solutions 
of the inverse band structure problem: for a given en- 
ergy E and k|| = they satisfy the Schrodinger equa- 
tion Hipk ± = Eipk ± in the bulk of the crystal. The cal- 
culations are performed with the ELAPW-kp method, 
which reduces the equation to a matrix eigenvalue prob- 
lem, with k± being the eigenvalues^ To the left from 
the matching plane zm, see Fig. (5J the LEED function 
$ is a linear combination of several ipk ± (only the waves 
with lmk± not exceeding 1 a.u. -1 are included). 

In the surface region, between zm and Zy, the poten- 
tial is different from the bulk potential, and the partial 
waves representation is not valid. Here the function $ is 
expanded in terms of the eigenfunctions £ n of the slab, 
which contains the scattering region, see Ref. HU- The 
functions £„, thus, have already taken into account the 
scattering by the overlayers. For each energy E the three 
representations are matched at the two planes zm and zy 
to construct a smoothly continuous function that satis- 
fies (with certain accuracy) the equation HQ = EQ in 
the embedded region. (The Schrodinger equation is sat- 
isfied by construction both in the bulk and in the vacuum 
half-spaces.) 

V. CLEAN Ni(lll) 

The ability of the substrate to conduct current can be 
characterized by the current carried by individual Bloch 
waves ijj. For the clean Ni(lll) surface the complex band 
structure in the TL direction and the energy-momentum 
distribution of the current are shown in Fig. \&[a.). The 
current carried (absorbed) by an individual Bloch wave 
is shown by the thickness of the dispersion curves. Below 
23 eV and above 34 eV relative to the Fermi energy the 
behavior of the main conducting branch is close to free- 
clcctron-likc, i.e., the electron transmission is effected 
by a single bulk band with almost parabolic dispersion. 
Between 23 eV and 34 eV, there is a wide Bragg gap 
centered at 28 eV and a number of special points, where 
the conducting branch switches from one bulk band to 
another. They give rise to the sharp peaks A, B, and C 
in the elastic (i.e., Vi = 0) R(E) spectrum of Ni(lll). 
In the present experiment only the peaks A and B are 
visible (in agreement with earlier work, Ref. [2~H ). This 
is explained by a strong effect of the inelastic scattering 
on the peak C. Figure^d) compares the measured R(E) 
spectrum (specular reflectivity) for clean Ni(lll) to the 
calculations with realistic values of the optical potential. 
The theoretical curve is the ratio of the specularly re- 
flected current to the incident current, and the experi- 
ment is the arbitrarily scaled /(V) curve. At moderate 
values of V ~ 1 eV the peak C completely disappears, in 
agreement with the experiment. The theoretical curves 
in Fig. EJd) are shifted by 1.6 eV to higher energies in 
order to bring the peaks A and B to the measured posi- 
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tions. This discrepancy is the result of the quasi-particle 
self-energy being inaccurately treated in the present cal- 
culation: the LEED functions are calculated as Kohn- 
Sham solutions in the local density approximation. The 
self-energy effect is seen to stretch the unoccupied spec- 
trum by 1.6 eV over the interval of 30 cV. 



VI. NiO(lll) CRYSTAL AND FILMS 

To illustrate the effect of a thin overlayer on the elec- 
tron reflection from Ni(lll), we have presented in Figs. [5] 
and[6^c) calculations for an oxygen monolayer fully com- 
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FIG. 6. (color online) (a) Conducting branches of the com- 
plex band structure (thick lines) superimposed onto the bulk 
band structure in the LTL interval (thin lines), (b) Normal 
incidence R{E) spectrum of Ni(lll) for Vi = 0.05 eV. (c) 
Normal incidence R(E) spectrum of an oxygen overlayer on 
Ni(lll), see Fig. [S] (d) Measured R(E) normal incidence 
spectrum of Ni(lll) (full line) and calculations (dashed and 
dot-dashed) for two choices of the energy dependence of the 
optical potential V\{E). Inset shows the two functions Vi(E) 
with respective line styles. The bars in the inset show the 
locations of the peaks A and B. 



mensurate with the substrate (thereby strongly com- 
pressively strained in the lateral direction if single-layer 
NiO(lll) is used as reference). Figure [5] compares the 
LEED functions for a clean Ni(lll) surface and for an 
oxygen overlayer on Ni(lll) at a kinetic energy of 10 cV. 
Already at the second Ni layer the wave function is seen 
to acquire the character of the substrate Bloch wave, but 
the oxygen overlayer reduces its amplitude. The compar- 
ison of the spectra in Figs. M^>) and [S{c) shows that in 
the interval 7-22 eV, where the band structure of the sub- 
strate is free-electron-like, the oxygen overlayer increases 
the specular reflectivity to as much as 50%. However, 
its presence does not result in any sharp structures in 
this energy interval: the reflectivity steadily increases at 
low kinetic energies similar to the case of clean Ni(lll). 
On the contrary, the measured spectrum of the oxidized 
surface, Fig. [TJd), shows a sharp decrease at low kinetic 
energies. 

Naturally, at low energies the electron reflection is 
very sensitive to the shape of the potential barrier be- 
tween the crystal half-space and the vacuum. This 
can be illustrated by a simple one-dimensional system, 
Fig. [TJa), in which the potential barrier is modeled by 
a two-step function. For a step-like potential barrier of 
height Vq [with V\ = 0, in the notation of Fig. [7|a)], the 
electron reflection R(E) steadily decreases with energy: 
R{E) = (V-E + V Q - VE~) 2 /(y/E + V a + s/E) 2 . However, 
when the width of the surface layer (with a different po- 
tential V\) is increased up to several atomic units, the 
R(E) curve becomes non-monotonic, and at low ener- 
gies transmission resonances appear - R(E) minima that 
characterize the surface layer. In the following, we will 
present an interpretation of the experimentally observed 
low-energy structure in the I(V) curve for an oxidized 
Ni(lll) surface as caused by a thin oxide layer. 
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FIG. 7. (color online) (a) Modeling of the overlayer by a 
two-step potential barrier. The constant bulk potential is 
Vo = 1 Ry. Effect of the surface layer potential Vi on the 
electron reflection for the width of the surface layer L = 2 a.u. 
(b), 5 a.u (c), and 9 a.u. (d). 
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FIG. 8. (color online) NiO trial structures on a model Ni(lll) 
substrate with the same Ni lattice as in NiO. From left to right 
one, two, and three layers of NiO(lll) on Ni(lll) as well as 
a truncated-bulk NiO(lll) crystal are shown. 
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FIG. 9. (color online) Measured R(E) normal incidence 
spectrum of 0/Ni(lll) (circles) and calculations for three 
(dashed) and two (solid line) surface layers of oxygen on a 
model Ni(lll) substrate with the same Ni lattice as in NiO. 
Theoretical curves are the ratio of the specularly reflected 
current to the incident current. They are shifted by 1.7 eV to 
higher energies to bring the main structures in agreement with 
their measured locations. The experiment is the arbitrarily 
scaled I(V) curve. The inset compares theoretical results for 
one, two, three layers of NiO, and for a NiO(lll) crystal. 



We have performed calculations of R(E) for one, two 
and three NiO layers on the model Ni substrate, as well 
as for bulk NiO(lll) (cf. Fig. [5]), which are shown in the 
inset of Fig. |HJ apart from the region of very low kinetic 
energies, the gross features of the R(E) curve are de- 
termined by the outermost NiO layers: the raise of the 
reflectivity below 12 eV is caused by the NiO overlayer, 
and the reflection from the Ni-NiO interface brings about 
only relatively weak wiggles. However, for the few-layer 
thin surface oxides, in contrast to bulk NiO, in approach- 
ing the vacuum level the reflectivity rapidly falls down to 
produce a very sharp minimum, in agreement with our 
measurements, where it is found at about 2.7 eV above 
the vacuum level. This dip is a signature of the con- 
ducting substrate: the bulk energy spectrum of NiO has 
a wide gap in this region, so the reflection by the bulk 
NiO crystal is complete, see the upper curve in the in- 
set of Fig. [5] (inelastic scattering is negligible at such low 
energies). The most interesting observation is that the 
presence of the conducting substrate does not lead to a 
uniform reduction of the reflectivity over the entire NiO 
band gap, but it leads to a sharp structure at a low kinetic 
energy. This minimum becomes deeper with decreased 
thickness of the NiO film, which resembles the transmis- 
sion resonances of the simple two-step model (Fig. [7j). 

The calculated specular reflectivity spectrum for two 
NiO layers agrees rather well with the measurements re- 
garding the overall shape, as well as the relative energy 
locations of main spectral structures: the discrepancies 
are within the expected error range due to possible errors 
in the quasiparticle self-energy The structures in the cal- 
culated spectrum are much sharper than in the measured 
curve, which may be attributed to the imperfect quality 
of the surface as well as to a finite angular and energy res- 
olution of the experiment. Additional discrepancies may 
also arise from the nanoscale heterogeneity of the system 
(cf. Sec. IIII B]) . which experimentally may lead to partial 
averaging over surface areas containing regions with only 
a single NiO layer. 



Because the formation of more than one NiO layer in 
the Ni(lll) lattice is improbable, we now consider the 
more realistic situation where the oxidized surface lay- 
ers have the structure of bulk NiO (see Fig. In view 
of the large lattice mismatch between the Ni substrate 
and the surface NiO layers, to construct a fully realis- 
tic model of the oxidized surface is beyond our compu- 
tational capabilities. However, we can make use of the 
observation that in the 15 eV wide energy region above 
the surface barrier, Ni(lll) has a free-electron-like band 
structure. Hence, we can model the underlying crystal by 
an expanded lattice of Ni, with the same lattice constant 
as NiO. Because the expanded Ni(lll) crystal too has a 
free-electron-like conduction band below E — Ep = 18 eV, 
this model appears quite plausible at low energies and 
capable of providing a qualitative understanding of the 
effect of the Ni substrate. 



VII. CONCLUSION 

We have studied the oxidation of the Ni(lll) surface 
in a concerted experimental and theoretical approach us- 
ing in- situ low-energy electron microscopy and ab ini- 
tio scattering theory. Upon exposure to molecular oxy- 
gen under ultrahigh- vacuum conditions, we observed the 
gradual formation of ultrathin NiO(lll) nuclei, which, 
at room temperature, form a quasi-continuous film. The 
resulting surface oxide exhibits a characteristic energy- 
dependent reflectivity R(E), which has been shown to 
provide detailed information on its electronic and geo- 
metric structure. 

For oxidation at room temperature, our results are ba- 
sically in line with previous publications by other groups 
reporting the formation of ultrathin surface oxides ex- 
hibiting thicknesses of two to three NiO layersj 7 ' 8 ' 10 At 
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elevated temperatures of 750 K, we have identified the 
growth of NiO(lll)-like patches, contrary to previous 
studies that have not suggested the growth of a locally 
wcll-dchncd surface oxide under these conditions. How- 
ever, in the present article we have also demonstrated 
that the evolving oxide "film" is highly heterogeneous 
with unreacted, essentially oxygen-free Ni(lll) areas in- 
between. An intriguing result of the paper is that a (lxl) 
periodicity is found for oxidation at high temperature 
where hydroxyl species should not be stable. Since we 
know from LEED investigations that the grain size is 
considerably larger for oxidation at 700-750 K instead of 
300 K, we conclude that neither a very small grain size 
nor termination by hydroxyls are necessary conditions for 
the stabilization of the unreconstructed (lxl) structure. 
While the (lxl) structure found at 300 K might still be 
induced by hydroxyls, this is not possible for oxidation 
at 750 K, leading us to the conclusion that its stability 
should rather be related to the finding of (i) an oxide sur- 
face unit cell that is rotated by 30° as compared to the 
Ni(lll) substrate, providing a changed interface to the 
underlying Ni and (ii) the ultrathin nature of the oxide 
film. 

Strong support for our findings is provided from the 
present ab initio scattering theory, which explains the 
observed normal incidence LEED spectra at very low en- 
ergies of both the clean and the oxidized Ni(lll) surface. 
In the latter case, the increase of the reflectivity in go- 
ing to lower energies is caused by a forbidden energy gap 
in the spectrum of the bulk NiO crystal along TL. The 
rapid decrease of R(E) at still lower energies appears 
only for a sufficiently thin NiO layer. Here, the charac- 
teristic feature is that the electron transmission through 
the oxide film does not increase uniformly over the gap 
with decreasing thickness of the oxide layer. Instead, a 
narrow transmission channel opens at the bottom of the 



gap, and it disappears again at the thickness of a single 
atomic NiO layer. These results indicate that in the ex- 
periment the oxide layer thickness is mostly two atomic 
layers, and that it can hardly be larger than three atomic 
layers. 

In essence, our experimental and theoretical findings 
support the formation of an ordered, ultrathin nickel 
oxide layer whose thickness is limited to less than one 
nanometer even at very high temperatures, while, in ac- 
cordance with ellipsometry data^ additional oxygen is 
randomly incorporated into the near-surface region upon 
prolonged oxygen exposure. This result implies that in 
the case of epitaxial oxide growth the general prediction 
by Cabrera and Mott may only hold for certain crys- 
tallographic orientations of the oxide film or the oxygen 
concentration profile in general, but does not necessarily 
allow conclusions on the actual thickness of the ordered 
oxide film. 
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